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ABSTRACT 


A  project  was  initiated  in  the  Spring  of  1960  to  calibrate  the 
Space  Surveillance  System  by  means  of  ballistic  cameras.  The 
principle  behind  this  calibration  is  a  comparison  between  the  posi¬ 
tions  of  the  Echo  balloon  as  determined  by  radio-frequency  reflec¬ 
tions  of  the  Space  Surveillance  System  and  as  determined  by  photo¬ 
graphs  against  a  star  background.  The  preliminary  results  show 
that  all  stations  have  zero  errors  of  1  e  s  s  than  0. 1  degree  at  the 
zenith.  A  second  error  having  a  standard  deviation  of  about  0.05 
degree  is  due  in  part  to  the  limited  precision  of  measurement  used. 
The  report  describes  the  work  done  prior  to  the  second  phase  of 
calibration,  that  phase  leading  to  a  much  higher  precision. 


PROBLEM  STATUS 

This  is  an  interim  report  on  the  first  phase  of  calibration; 
work  is  continuing. 
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OPTICAL  CALIBRATION  OF  THE  U.S.  NAVAL  SPACE 
SURVEILLANCE  SYSTEM 


INTRODUCTION 

The  U.S.  Naval  Space  Surveillance  System  contains  four  radio  interferometer-type 
receiving  stations  located  along  a  great  circle  which  crosses  the  southern  part  of  the 
United  States.  The  receiving  stations  are  located  at  Ft.  Stewart,  Georgia;  Silver  Lake, 
Mississippi;  Elephant  Butte,  New  Mexico;  and  Brown  Field,  Chula  Vista,  California.  Three 
transmitters  are  located  at  suitable  points  remote  from  the  receiving  stations  along  the 
same  great  circle.  These  transmitters  provide  sources  of  electromagnetic  radiation 
which  induce  secondary  radiation  from  satellites  when  they  enter  the  fields  generated 
by  the  transmitting  stations. 

The  transmitting  and  receiving  antennas  form  fan- shaped  patterns  along  the  plane  of 
the  great  circle,  so  that  the  region  of  detection  is  confined  approximately  to  this  plane. 
Hence,  an  angular  measurement  of  the  direction  to  a  satellite  taken  simultaneously  from 
each  of  two  receiving  stations  is  sufficient  to  determine  the  three  position  coordinates  of 
the  satellite  at  the  time  of  measurement.  The  purpose  of  optical  calibration  is  to  deter¬ 
mine  any  small  errors  existing  in  the  angular  measurements  obtained  with  the  inter¬ 
ferometer,  and  thus  improve  the  accuracy  of  the  system.  Since  the  far  field  of  the 
antennas  is  higher  than  conventional  aircraft  fly,  the  method  chosen  uses  photographs  of 
the  Echo  balloon  against  a  star  background.  Simultaneous  radio-frequency  reflection 
from  the  satellite  permits  the  system  to  be  calibrated. 

In  anticipation  of  the  launching  of  the  Echo  satellite,  a  project  to  calibrate  the  Space 
Surveillance  System  by  optical  means  was  assigned  in  April  1960.  A  minimum  accuracy 
of  measurement  of  0.1  degree  was  considered  sufficient  to  provide  useful  calibration  data 
and  was  adopted  as  an  initial  goal.  To  carry  out  the  project  it  was  necessary  to  acquire 
and  install  suitable  cameras  at  each  receiving  station,  and  to  provide  these  . cameras  with 
time  coding  equipment.  Also,  a  reduction  procedure  was  necessary  to  obtain  the  position 
of  the  satellite  with  respect  to  the  receiving  station  as  a  function  of  time,  as  the  angle  so 
determined  optically  would  be  compared  with  the  angle  measurement  taken  from  the 
interferometer  to  obtain  the  calibration  error. 


INSTRUMENTATION 

At  the  time  this  project  was  initiated  a  camera  employing  a  “K38”  f/6.3  lens  with  a 
24-inch  focal  length  and  an  equatorial  mount  with  sidereal  drive  was  available  at  the  NRL 
Chesapeake  Bay  Annex.  Time  coding  equipment  devised  by  E.  Habib  of  NASA  was  used 
with  this  camera.  Habib’s  method  involved  displacing  the  photographic  plate  with  a  sole¬ 
noid  actuated  plunger.  After  making  initial  tests  on  this  camera,  the  entire  assembly  was 
shipped  to  the  Ft.  Stewart,  Georgia,  Space  Surveillance  Station  and  put  into  operation  in 
October  1960.  Figure  1  shows  a  photograph  taken  with  this  camera. 

The  breaks  in  the  satellite  trail  produced  by  the  time  coding  are  clearly  visible  on 
the  trail  of  the  satellite.  The  time  code  marks  appear  at  2-secohd  intervals  except  during 
the  readout  of  Universal  Time,  which  starts  exactly  on  each  minute.  The  digital  time  code 
is  formed  by  pulses  separated  by  1/5  second.  Four  bursts  of  pulses  appear  at  2-second 
intervals  to  give  the  four  digits  required  to  read  the  time  of  day  on  a  24-hour  clock.  As 
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a  consequence  of  the  plate  displacement  during  exposure,  the  stars  have  double  images. 

The  displaced  images  are  much  fainter  because  of  the  relatively  small  exposure  time  and 
cause  no  problem  in  plate  reduction. 

A  small  plywood  box  camera  loaned  by  Norton  Goodwin  of  the  Volunteer  Satellite 
Tracking  Program  was  set  up  at  Silver  Lake,  Mississippi.  It  employed  an  f/2.5  Aero 
Ektar  lens  with  a  7-inch  focsil  length.  Despite  apparent  crude  construction  of  this  camera, 
its  optical  quality  was  excellent.  This  camera  was  used  to  obtain  trail  type  images. 

Breaking  of  the  satellite  trail  was  done  manually  with  a  mask,  a  method  used  to  some 
extent  by  amateur  satellite  photographers.  This  camera  provided  useful  data  as  early  as 
September  12,  1960,  at  the  Silver  Lake  site. 

Construction  of  three  additional  equatorial  mounts  of  the  type  used  at  Ft.  Stewart, 
Georgia  was  begun  in  the  Fall  of  1960.  The  new  mounts  were  equipped  with  K37-type 
aircraft  cameras.  The  shutter  actuation  mechanism  in  these  cameras  was  so  designed 
that  the  shutters  could  be  opened  by  the  application  of  28  volts  dc.  To  utilize  this  built-in 
capability  a  notched  disc  was  fabricated  and  mounted  on  a  1-rpm  synchronous  motor  shaft. 
The  notches  alternately  opened  and  closed  a  microswitch  in  the  28-volt  line  to  the  shutter 
drive  mechanism.  This  voltage  was  also  fed  to  a  multichannel  recorder  and  displayed  with 
Universal  Time  and  other  information.  Identifying  notches  were  inserted  at  120-degree 
intervals  around  the  notched  disc,  making  each  20-second  interval  unique  within  a  1-minute 
interval.  Since  auxiliary  information  taken  from  the  interferometer  records  defines  the 
time  at  which  the  satellite  crossed  the  antenna  beam,  no  problem  of  identifying  the  minute 
interval  exists.  Figure  2  is  an  example  of  the  type  of  time  coding  obtained  by  this  method. 

Figure  3  is  a  sample  of  recorded  data  which  gives  the  phase  information  from  which 
the  zenith  angle  of  the  satellite  is  determined.  Along  the  lower  margins  appears  the  coded 
Universal  Time.  Directly  above  the  time  the  age  response  is  shown  as  the  satellite  crossed 
the  plane  of  the  great  circle  usually  referred  to  as  the  “fence.”  In  the  next  channel  is  a 
record  of  the  voltage  applied  to  the  shutter  drive  mechanism.  This  is  used  in  determining 
the  Universal  Time  for  points  along  a  satellite  trail  such  as  the  one  shown  in  Fig.  2.  The 
installation  of  the  three  additional  cameras  at  the  receiving  sites  was  completed  in  May  1961. 


DATA  REDUCTION 

The  modest  requirement  of  0.1-degree  accuracy  in  angle  determination  allowed  some 
choice  of  method  in  reducing  photographic  plates.  The  method  used  is  essentially  as 
follows:  A  star  chart  photographed  to  the  same  scale  as  the  plates  to  be  reduced  is  placed 
over  an  illuminated  viewer.  The  photographic  plate  to  be  reduced  is  placed  over  the  star 
chart  and  adjusted  to  match  the  stars  on  the  chart.  The  path  of  the  satellite  trail  among 
the  stars  is  carefully  observed  and  plotted  on  the  chart.  Time  information  may  also  be 
designated  along  the  satellite  trail  by  reference  to  the  photographic  plate.  Figure  4  is  a 
plot  of  such  a  satellite  trail.  It  is  also  necessary  to  find  the  projection  among  the  stars 
of  the  fence.  Since  the  celestial  coordinates  are  based  on  the  position  of  the  vernal  equinox, 
there  is  a  slow  but  continual  drift  of  the  coordinate  grid  among  the  stars.  To  compensate 
for  this  the  present-day  celestial  coordinates  of  the  fence  are  transformed  to  the  set  used 
in  preparing  the  star  chart.  (In  the  case  of  the  Bonner  Durchmusterung  the  Epoch  date  for 
the  chart  coordinates  is  1855.)  The  information  on  the  fence  coordinates  is  plotted  in  the 
neighborhood  of  the  satellite  trail  which  intersects  it.  One  degree  intervals  are  marked 
along  the  line  representing  the  fence  in  Fig.  4.  The  displacement  of  the  time  of  fence 
crossing  indicated  by  the  interferometer  system  from  that  obtained  by  the  optical  system 
has  been  too  small  to  be  considered  in  the  method  described  above.  On  the  other  hand  the 
differences  in  the  optically  measured  zenith  angles  and  the  angles  obtained  from  the  inter¬ 
ferometer  have  been  quite  appreciable.  It  is  this  zenith  angle  error  which  has  been  inves¬ 
tigated  in  some  detail.  Measurements  of  optical  zenith  angle  are  read  to  0.01  degree  from 
the  chart.  This  graphical  procedure  has  little  of  the  elegance  of  the  method  using  an  optical 


4 


NAVAL  RESEARCH  LABORATORY 


Fig.  2  -  Example  of  an  Echo  satellite  trail  showing  time  coding 
produced  by  an  electrically  actuated  shutter  mechanism 
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SHUTTER  DRIVE  VOLTAGE  . 


Fig.  3  -  Sample  of  the  recorded  data  from  which  the  zenith  angle 
of  the  satellite  is  determined 
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comparator;  but  neither  does  it  require  as  much  time  to  obtain  reasonably  good  informa¬ 
tion.  It  has  been  estimated  that  the  overall  zenith  angle  measurement  error  is  about  0.05 
degree.  This  error  is  attributed  largely  to  plotting  and  reading  of  the  graph.  The  data 
reduced  by  this  method,  which  is  treated  below,  verifies  that  this  is  a  reasonable  estimate 
of  the  error  of  this  method. 


ANALYTICAL  METHODS 

It  has  been  indicated  above  that  only  the  zenith  angle  error  is  of  consequence.  The 
remainder  of  this  report  is  concerned  with  comparing  the  zenith  angle  obtained  by  optical 
methods  with  angle  measurements  taken  from  the  interferometers  and  the  data  transmis¬ 
sion  system.  Figure  5  presents  a  block  diagram  of  one  pair  of  antennas  used  in  the  radio 
interferometer  and  traces  the  signal  from  the  receiving  site,  to  the  data  processing  center 
at  the  Naval  Weapons  Laboratory,  Dahlgren,  Virginia,  where  the  phase  information  from 
the  antenna  pair  appears  on  one  channel  of  a  multichannel  recording  similar  to  the  one 
shown  in  Fig.  3.  The  phase  data  recorded  at  NWL  has  been  reduced  at  NRL,  using  in  part, 
calibration  data  readings  made  at  NWL.  A  summary  of  the  reduction  of  these  data,  is 
given  in  Tables  1-4  at  the  end  of  the  report.  (The  details  of  these  tables  will  be  discussed 
later  in  the  report.)  A  second  source  of  data  has  been  the  records  taken  at  the  receiving 
sites,  which  is  in  part  a  duplication  of  the  NWL  data. 
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Fig.  5  -  The  processing  of  the  signal 
from  one  pair  of  antennas 

By  comparing  the  phase  data  taken  at  these  two  points  in  the  system,  it  has  been 
possible  to  look  for  any  degradation  in  the  phase  information  during  its  transmission  to 
the  remote  data  center  at  NWL.  The  two  sources  of  recorded  data  also  provided  a  check 
on  time  coding  and  possible  errors  in  the  reduction.  The  angle  information  obtained  from 
the  two  sources  has  been  compared  with  optical  measurements. 

In  making  this  comparison,  it  was  found  desirable  to  depart  slightly  from  the  cus¬ 
tomary  angle  designation  used  in  Space  Surveillance.  Customarily  the  direction  of  arrival 
of  a  signal  from  a  satellite  is  described  as  6  degrees  east  or  6  degreeswestof  zenith  accord¬ 
ing  to  the  case.  This  is  equivalent  to  using  spherical  coordinates  in  which  the  zenith  angle  e 
plays  the  conventional  role  and  4.  takes  one  of  two  values  separated  by  180  degrees.  It  was 
much  more  convenient  for  this  analysis  to  use  a  single  coordinate  as  shown  in  Fig.  6. 


O 

ZENITH 


Fig.  6  -  Single -angle  coordinate  ranging 
from  -90  to  +90  degrees 
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This  was  done  by  assigning  a  positive  value  to  zenith  angles  in  the  east  and  a  negative 
value  to  those  in  the  west.*  By  doing  this  d  as  well  as  sin  0  are  continuous  single-valued 
functions  throughout  the  region  of  interest. 

Perhaps  the  most  natural  form  in  which  to  present  errors  in  the  interferometer  angle 
measurement  is  in  terms  of  angle  error  as  a  function  of  zenith  angle.  This  has  been  done 
by  plotting  the  difference  function 


for  Ft.  Stewart,  Georgia,  in  Fig.  7(a)  and  for  Silver  Lake,  Mississippi,  in  Fig.  7(b),  with 
6,  being  defined  as  the  angle  measured  by  the  interferometer  and  e  that  determined  by 
optical  instruments.  ° 

In  Fig.  7(a)  it  is  seen  that  a  marked  increase  in  angular  error  takes  place  for  zenith 
angles  beyond  -70  degrees  in  the  west.  A  curve  having  the  best  fit  passed  through  the 
distribution  of  points  may  be  used  to  obtain  a  mean  value  of  error  in  the  interferometer 
data  as  a  function  of  zenith  angle.  This  procedure,  although  mathematically  valid  provides 
no  insight  into  the  cause  of  the  error  and  fails  to  provide  a  firm  basis  for  projecting  the 
curve  over  angles  for  which  no  optical  observations  have  been  made.  Figure  7(b)  presents 
a  more  severe  condition  which  would  be  difficult  to  resolve  without  resorting  to  a  more 
powerful  method  of  analysis.  To  continue  the  analysis,  one  seeks  an  analytical  expression 
for  a  curve  that  will  fit  the  observed  points  of  Fig.  7(a),  or,  what  is  equivalent,  a  trans¬ 
formation  in  which  the  curve  through  the  points  degenerates  to  a  straight  line. 

In  Fig.  8(a),  the  data  in  Fig.  7(a)  has  been  plotted  in  the  form 


Asin^  -  sin^,  -  sinfi  . 

*  o 


(2) 


It  is  apparent  that  the  desired  straight  line  has  been  found.  The  resulting  simple  form 
permits  data  to  be  summarized  by  statistical  averaging.  This  has  been  done  for  all  base¬ 
lines  and  the  results  are  given  in  Tables  1-4.  The  error  in  the  sine  of  the  zenith  angle  is 
independent  of  the  zenith  angle  within  the  limits  of  the  accuracy  of  the  equipment  used  for 
making  optical  measurements.  Since  the  dependence  on  angle  of  arrival  was  removed  by 
the  transformation  to  sine  functions,  it  is  useful  to  replace  the  zenith  angle  by  time  as  a 
plotting  parameter,  as  the  possibility  of  calibration  changes  over  long  periods  of  time  is 
one  of  the  questions  for  which  an  answer  has  been  sought.  Accordingly  the  plot  for  Silver 
Lake  is  presented  as  a  function  of  time  in  Fig.  8(b). 

The  time  dependent  plot  reduces  the  confusion  of  Fig.  7(b)  to  order.  It  is  apparent  that 
a  new  factor  entered  the  measurements  at  Silver  Lake.  A  change  in  the  calibration  error 
appears  to  have  taken  place  between  January  10  and  12,  1961.  (A  second  shift  in  calibra¬ 
tion  was  introduced  between  1800  and  2020Z  June  28,  1961,  to  correct  the  condition  which 
existed  during  the  previous  months.)  To  explain  the  nature  of  the  error  which  has  been 
found  in  the  graphs,  let  us  refer  to  the  basic  geometry  of  the  system.  In  Fig.  9  is  shown 
a  pair  of  antenna  arrays  Aj  and  Aj  separated  by  a  distance  B.  6  is  the  angle  between  the 
zenith  and  the  direction  of  arrival  of  a  signal.  It  is  clear  that 

(3) 


♦The  choice  of  signs  used  here  was  suggested  by  an  earlier  unpublished  paper  by  the  author 
in  which  dwasdefined  as  the  angle  between  the  easternhorizon  and  the  direction  of  arrival  of 
the  signal.  The  range  of  6  was  0  <  Q  <180°,  and  cos  6  in  that  treatment  had  the  same  value 
as  sin  6  does  here. 
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Fig.  7(a)  -  Zenith  angle  error.  Ft.  Stewart,  Georgia 


Fig.  7(b)  -  Zenith  angle  error.  Silver  Lake,  Mississippi 
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Fig.  8(a)  -  Data  of  Fig.  7(a)  plotted  in  the  form 
of  Eq.  (2)  as  a  function  of  zenith  angle 
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Fig.  8(b)  -  Data  of  Fig.  7(b)  plotted  in  the  form 
of  Eq.  (2)  as  a  function  of  time 
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where  n  is  an  integer  and  4>  (0  <  4>  <  l)  is  the  fractional  part  of  a  wavelength.  Since  the  phase 
detector  n  measures  the  difference  in  phase  at  its  inputs  and  a^  rather  than  the  phase 
equivalent  of  the  distance  s,  any  difference  in  the  phase  length  0^  of  the  transmission 
lines  to  Aj  and  Aj  appears  in  the  output  as  an  additive  constant  to  n  +  0.  This  error 
converts  the  theoretical  Eq.  (3)  obtained  by  geometric  considerations  into 

sin<9  +  Asine  =  +  Is.  =  m  f  .  (4) 

A  plot  of  Eq.  (3)  (solid  lines)  is  given  in  Fig.  10.  It  is  seen  that  (n  +  0)/b  is  a  straight  line 
passing  through  the  origin  at  the  zenith  and  is  a  composite  of  the  stairstep  function  n  and 
the  sawtooth  function  0 .  A  second  plot  in  dashed  lines  has  been  made  for  (m  +  $)/3 
illustrating  the  linear  relationship  between  the  two  functions.  It  is  clear  from  Fig.  10 
and  Eq.  (4)  that  the  desired  phase  0  may  be  extracted  from  the  observed  phase  $, 
provided  one  knows  the  value  of  0,  .  From  Eqs.  (3)  and  (4)  one  obtains 

0g  =  B  Asin^  .  (5) 


5  — 

5  — 

4  — 

3  — 

Z  — 

- ^ 

\ 

1  1  ^1111 

1  1  1  1  1  1  1 _ L  1  - 

-1.0  -0.9  -0.8  -0.7  -0.6  -0,5  -04  -03  -0.2  -01  0.0  0.1  0.2  03  0.4  0.5  0.6  07  08  09  1.0 


-90*  -30*  SIN  e  +30*  -►90* 

WEST  EAST 


Fig.  10  -  Plots  of  Eqs.  (3)  and  (4) 
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In  practice  one  must  recognize  that  the  direction  of  arrival  e  is  not  measured  exactly  and 
that  the  phase  detector  may  introduce  a  random  phase  error.  To  obtain  a  more  precise 
value  for  one  must  take  the  mean  of  a  number  of  observations.  Figure  11  shows  a 
plot  of  the  convergence  of  the  mean  value  of  Asin0  as  a  function  of  the  number  of  measure¬ 
ments,  n  •  On  the  same  graph  the  function 


has  been  plotted,  a  is  obtained  from  the  complete  set  of  observations  plotted  on  the  graph. 

is  taken  as  an  adequate  approximation  for  the  probable  error  in  the  mean  and  is 
included  in  three  proportional  forms  in  Tables  1-4  for  maximum  convenience  for  different 
applications.  Because  of  the  small  values  of  probable  error  in  where  an  adequate 
number  of  observations  have  been  made,  it  is  apparent  that  this  error  can  be  eliminated 
effectively.  This  may  be  done  in  two  ways: 

1.  may  be  converted  to  inches  of  transmission  line  and  a  suitable  correction 
section  added  (or  subtracted)  according  to  the  case  at  a  convenient  place. 

2.  4>^  X  100  is  the  reading  conventionally  taken  from  the  recorders;  hence  x  100 
may  be  subtracted  from  the  recorder  output  readings  without  making  any  physical  change 
in  the  system. 

Method  1  has  the  advantage  that  once  made  it  requires  no  further  attention.  Since  a 
number  of  weeks  is  required  to  recalibrate  the  system  optically,  the  new  calibration  con¬ 
stant  cannot  be  checked  immediately. 

Method  2  is  applicable  to  past  data  and  can  be  used  to  study  the  improvement  that 
results  from  application  to  previous  analyses.  An  example  of  this  is  illu  trated  when  the 
corrections  given  in  Tables  1-4  are  applied  to  the  data  from  Ft.  Stewart,  Georgia,  and 
Silver  Lake,  Mississippi.  When  this  is  done  the  angle  errors  shown  in  Figs.  7(a)  and 
7(b)  reduce  to  those  shown  in  Figs.  12(a)  and  12(b).  In  this  case  the  improvement  brought 
about  by  applying  corrections  is  quite  impressive. 


Fig.  11  -  Variation  of  the  mean  error  in  sin  6  vs  the 
number  of  observations  (Ft.  Stewart) 
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Fig.  12(a)  -  Residual  errors  in  the  Ft.  Stewart  data 
of  Fig.  7(a)  after  applying  the  corrections  of  Table  1 


Fig.  12(b)  -  Residual  errors  in  the  Silver  Lake  data 
of  Fig.  7(b)  after  applying  the  corrections  of  Table  2 
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The  standard  deviation  a,  which  measures  the  probable  error  of  an  individual  obser¬ 
vation,  is  a  more  persistent  error  of  the  system.  However,  due  to  the  elaborateness  of 
the  system  a  partial  analysis  can  be  made  of  this  variability  which  leads  to  additional 
information  on  the  capability  of  the  system.  It  is  well  to  proceed  by  considering  some 
possible  causes  of  the  scatter  in  measurements,  which  are: 

1.  Random  shifts  in  the  phase  measuring  equipment  and  data  transmission  system. 

2.  Random  errors  in  reading  recording  data. 

3.  Random  departures  of  the  direction  of  arrival  of  radio  signals  from  the  optically 
observed  direction. 

4.  Random  errors  in  the  measurements  of  optical  angles. 

These  four  sources  of  variations  divide  into  two  groups  depending  on  whether  the 
error  is  in  the  measurement  of  4>  or  sin  e .  If  one  considers  the  equation 


i4>  -  B  Ssiii0  (7) 

one  sees  that  a  variation  in  s<^■  due  to  a  variation  in  Ssinfi  approaches  zero  linearly  as 
R  approaches  zero.  Hence  the  standard  deviation  of  &4>  is  independent  of  variability  in 
the  angle  error  Bsine  for  a  sufficiently  short  baseline.  Examination  of  the  values  in  the 
tables  for  i4>  shows  no  trend  for  baselines  not  exceeding  136  ft.  One  is  led  to  assume  that 
the  differences  in  the  h4>  values  for  the  various  baselines  results  from  a  choice  of 
samples  and  variability  in  quality  of  different  sets  of  electronics.  The  average  value  of 
the  standard  deviation  of  S<p  for  NWL  data  for  all  these  baselines  is  0.0394\.  Using  this 
value  one  may  derive  a  capability  limit  of  the  system  for  a  given  baseline  by  substituting 
this  value  into 


8sin5  =  (8) 

Solving  this  equation  indicates  that  a  520-ft  baseline  will  have  an  equivalent  zenith  angle 
error  of  0.04  degree  because  of  internal  system  variations.  Since  other  random  effects 
may  appear  as  the  precision  of  calibration  is  extended,  it  is  not  known  as  yet  how  much 
the  actual  random  errors  will  exceed  these  known  limitations  due  to  internal  system 
effects.  A  larger  standard  deviation  is  obtained  for  the  three  stations  using  the  520-ft 
baselines  than  is  obtained  from  the  shorter  baselines  at  the  same  stations.  Since  we  have 
no  reason  to  suspect  these  baselines  to  have  more  noise  than  the  shorter  ones,  this 
increase  must  be  attributed  to  a  variation  in  the  measured  value  of  e  due  to  causes  3  and 
4  listed  above.  Since  the  relationship  of  these  variations  is  completely  random  in  respect 
to  the  variation  (0.0394^)  for  phase  noise  &4>,  one  may  extract  a  value  for  the  variation  due 
to  variation  in  the  value  of  e  by 


Ssin0  -lASsinS)^  +  (SsinS)^  (9) 

f  9  <t> 

in  which  Ssine  is  the  standard  deviation  given  under  AsinS  in  the  tables,  ( Ssinfl)  is  the 
contribution  to  the  variability  in  8  sine  computed  for  the  noise  figure  (0.0394\)  in  the 
phase  measuring  equipment,  and  (8sine)g,  which  is  to  be  determined,  is  the  variability  in 
the  optical  measurement  of  6.  Solution  of  Eq.  (9)  gave  a  value  of  0.057  degree  for  the 
equivalent  zenith  angle  error  due  to  variations  in  e  .  Since  the  rise  in  standard  deviation 
was  not  evident  at  Ft.  Stewart,  where  a  higher  quality  optical  system  is  in  use,  it  is 
believed  that  this  error  can  be  attributed  largely  to  limited  capability  of  the  optical 
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measurements.  The  value  0.057  degree  derived  from  statistics  on  the  system  agrees 
quite  well  with  the  figure  0.05  degree  mentioned  earlier  which  was  made  as  an  intelligent 
guess  before  any  of  this  reduction  had  been  made. 


DISCUSSION  OF  GRAPHS  AND  TABLES 

Some  additional  details  concerning  the  graphs  and  tables  remain  to  be  discussed. 
Figures  7(a)  and  7(b)  give  angle  error  plots  for  the  two  stations  in  the  East  for  which  the 
largest  number  of  observations  have  been  made.  A  theoretical  curve  is  drawn  in  Fig.  7(a) 
which  is  a  transformation  of  the  straight  line  mean  given  in  Fig.  8(a).  The  agreement 
between  this  curve  and  observations  is  quite  good.  Similar  curves  are  drawn  in  Fig.  7(b) 
for  the  three  time-dependent  mean-value  lines  of  Fig.  8(b).  The  zenith  angle  plot  used  in 
Fig.  8(a)  has  been  replaced  by  a  time-dependent  plot  in  Fig.  8(b).  To  present  all  of  the 
data  for  the  period  on  one  graph,  gaps  in  the  time  scale  indicated  by  dark  dashes  have 
been  inserted.  An  equivalent  zenith  angle  error  of  0. 1  degree  is  also  plotted  in  Figs.  8(a) 
and  8(b).  The  equivalent  zenith  angle  is  defined  as  follows: 


Equiv.  zenith  angle  -  arcsin  |Asin0|  . 


The  residual  angle  errors  for  Ft.  Stewart  and  Silver  Lake  were  obtained  by  extracting  the 
proper  zero  or  bias  error  from  the  phase  readings  from  these  two  stations.  These  plots 
in  Figs.  12(a)  and  12(b)  are  the  appropriate  form  for  investigating  the  effects  of  refraction. 
It  appears  from  the  data  shown  here  that  the  measurement  of  refraction  is  beyond  the 
capabilities  of  the  data  reduction  method  employed.  Graphs  for  Elephant  Butte  and 
San  Diego  have  not  been  given  since  ali  information  of  interest  can  be  presented  more 
conveniently  in  tabular  form. 

The  tabulated  data  have  been  arranged  to  give  convenient  values  for  various  types  of 
applications.  AslniS (which  is  the  mean  of  n  observations)  and  its  standard  deviation  were 
the  first  values  extracted  by  the  analytical  process  employed.  In  Figs.  8(a)  and  8(b)  the 
appropriate  vaiues  of  Asini?  have  been  inserted  and  identified  as  the  mean  value  for  the 
distributions  of  points.  In  the  tables  the  value  cr/>/n  is  also  given  as  a  measure  of  the 
probable  error  in  the  mean.  Asing  is  converted  to  its  equivalent  in  terms  of  by  Eq.(5). 
This  is  given  in  two  forms:  the  first  in  terms  of  degrees  is  convenient  for  comparison 
with  electrical  measurements  of  line  length  give  results  directly  in  electrical  degrees; 
the  second  form  of  4‘^  is  equai  to  \  x  100  so  that  the  error  so  quoted  is  equal  to  percentage 
error  in  defiection  on  the  data  recorders.  The  corresponding  values  cr//n  of  precision  in 
are  given  as  a  confidence  factor  for  the  value  of  4>^.  The  correction  factors  for  base¬ 
lines  other  than  the  largest  have  no  effect  on  the  final  value  of  e,  provided  they  do  not 
cause  an  error  in  the  value  of  n  of  Eq.  (3).  However  they  may  have  a  use  in  connection 
with  data  analysis  and  in  setting  up  automatic  angle  resolving  equipment  and  accordingly 
have  been  included.  In  some  cases  definite  shifts  in  calibration  have  been  observed  and 
these  have  been  identified  in  the  tables. 


SUMMARY 

A  preliminary  calibration  of  the  Space  Surveillance  System  has  met  the  initial  require¬ 
ment  for  a  precision  of  0. 1  degree.  A  bias  caused  by  an  error  in  effective  transmission 
line  length  has  been  identified.  This  can  be  eliminated  by  either  a  mathematical  or  physi¬ 
cal  correction  to  the  system.  A  second  type  of  error  is  measured  in  terms  of  standard 
deviation.  A  component  of  this  error  has  been  traced  to  random  internal  phase  shifts  in 
the  interferometer  system.  A  second  component  due  to  random  variations  in  the  apparent 
angle  of  arrival  are  caused  in  part  by  limited  precision  in  the  optical  reduction  procedure. 
Comparison  of  data  taken  from  records  at  the  receiving  sites  and  NWL  show  no  significant 
changes  in  the  data  transmission  process. 
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More  refined  data  reduction  procedures  are  being  initiated  to  reduce  the  errors 
attributed  to  the  measurements  of  the  optical  positions. 
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Table  1 

Ft.  Stewart,  Georgia,  Calibration 


Baseline 

Length 

(ft) 

No. 

of 

Obs. 

Error  Expressed  as  a 
Function  of  sine 

Error  Expressed 

as  d’e  In 
Terms  of  Elec¬ 
trical  Degrees 

Error  Expressed 
as  in  Terms  of 
Wavelength* 

(%) 

Equiv. 

Zenith 

Angle 

Error 

(deg) 

Mean  of 

Asin6^ 

Std.  Dev., 

O' 

(  Ssin6^) 

Probable 
Error  in 
the  Mean, 

cr//n 

Mean 

Probable 
Error  in 
Mean, 

cr//n 

Mean 

Probable 
Error  in 
Mean, 

a//n 

413 

51 

-0.0014 

±0.0009 

.....  — 

±0.0001 

-22.4 

±2.2 

-  6.2 

±0.6 

00 

o 

o 

1 

±0.05 

136 

51 

-0.0067 

±0.0027 

±0.0004 

-36.2 

±2.0 

-10.1 

52 

51 

-0.0021 

±0.0049 

±0.0007 

-  4.3 

±1.4 

■h 

±0.4 

28 

50 

-0.0032 

±0.0083 

±0.0012 

-  3.5 

±1.3 

±0.3 

24 

50 

+0.0061 

±0.0154 

±0.0022 

+  5.8 

±2.1 

■g 

±0.6 

16 

51 

-0.0174 

±0.0174 

±0.0025 

-11.0 

±1.5 

BB 

±0.4 

*100%  recorder  deflection  is  equal  to  one  wavelength  in  terms  of  phase;  hence  the  numbers  in  this 
column  are  equal  to  the  percentage  of  recorder  deflection. 


Table  2 

Silver  Lake,  Mississippi,  Calibration 


Baseline 

Length 

(ft) 

No. 

of 

Obs. 

Error  Expressed  as  a 
Function  of  sin6 

_ 

Error  Expressed 
as  0,  in 
Terms  of  Elec¬ 
trical  Degrees 

Error  Expressed 
as  d>,in  Terms  of 
Wavelength 
(%) 

Equiv. 

Zenith 

Angle 

Error 

(deg) 

Mean  of 

Asin5 

Std.  Dev., 

cr 

(  Ssin^') 

Probable 
Error  In 
the  Mean, 

u/'/n 

Mean 

Probable 
Error  in 
Mean, 

<j//n 

Mean 

Probable 
Error  in 
Mean, 

a/Zn 

521* 

9 

-0.0007 

±0.0015 

±0.0005 

-15.3 

±10.1 

-4.2 

±2.8 

±0.085 

129t 

19 

+0.0017 

±0.0027 

±0.0006 

+6.6 

±3.2 

+2.4 

±0.9 

60.8 

27 

+  0.0024 

±0.0168 

±0.0032 

+5.9 

±7.8 

+  1.6 

±2.2 

52.8 

56 

+0.0014 

±0.0093 

±0.0012 

+3.2 

±2.6 

+0.8 

±0.7 

20 

56 

+0.0025 

±0.0169 

±0.0023 

±1.8 

+0.6 

±0.5 

16 

56’ 

-0.0111 

±0.0203 

±0.0027 

-7.0 

±1.7 

-1.9 

±0.5 

5211 

33 

+0.0061 

±0.0011 

±0.0002 

+125.2 

±4.1 

+34.8 

±1.1 

±0.065 

129§ 

37 

-0.0013 

±0.0017 

±0.0003 

-6.5 

±1.3 

-1.8 

±0.4 

520T 

12 

-0.0017 

±0.0015 

±0.0004 

-35.1 

±8.9 

-9.7 

±2.5 

boHiI 

±0.085 

^Applies  from  Sept.  12,  I960,  to  Jan.  10,  1961. 
tApplies  from  Sept.  12,  I960,  to  Feb.  15,  1961. 
TApplies  from  Jan.  12,  1961,  to  June  28,  1961. 
§Applies  from  Mar.  9,  1961,  to  Oct.  21,  1961. 
TApplies  after  June  28,  1961. 
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Table  3 


Elephant  Butte,  New  Mexico,  Calibration 


Baseline 

Length 

(ft) 

No. 

of 

Obs. 

Error  Expressed  as  a 
Function  of  sin^ 

Error  Expressed 
as  in 
Terms  of  Elec¬ 
trical  Degrees 

Error  Expressed 
as  <Pein  Terms  of 
Wavelength* 

(%) 

Equiv. 

Zenith 

Angle 

Error 

(deg) 

Mean  of 

Asin^ 

Std.  Dev., 

a 

(Ssin^?) 

Probable 
Error  in 
the  Mean, 

<r//n 

Mean 

Probable 
Error  in 
Me  am 

a//n 

Mean 

Probable 
Error  in 
Meam 

o-//n 

520 

-0.0001 

±0.0007 

±0.0002 

-1.0 

±4.8 

-0.3 

±1.3 

0.00 

±0.04 

128 

+  0.0021 

±0.0026 

±0.0008 

+10.7 

±4.1 

+3.0 

±1.1 

60 

+0.0015 

±0.0066 

±0.0021 

+3.6 

±4.9 

+1.0 

±1.4 

52 

+0.0041 

±0.0088 

±0.0028 

+8.4 

±5.8 

+2.3 

±1.6 

20 

9 

+0.0023 

±0.0086 

±0.0029 

+1.8 

±2.3 

+0.5 

16 

+0.0178 

±0.0232 

±0.0073 

+  11.3 

±4.6 

+3.1 

±1.3 

“i'100%  recorder  deflection  is  equal  to  one  wavelength  in  terms  of  phase;  hence  the  numbers  in  this 
column  are  equal  to  the  percentage  of  recorder  deflection. 


Table  4 

San  Diego,  California,  Calibration 


Baseline 

Length 

(ft) 

No. 

of 

Obs. 

Error  Expressed  as  a 
Function  of  sin0 

Error  Expressed 
as  4:^  in 
Terms  of  Elec¬ 
trical  Degrees 

Error  Expressed 
as  in  Terms  of 
Wavelength 
(%) 

Equiv. 

Zenith 

Angle 

Error 

(deg) 

Mean  of 
Asin(? 

Std.  Dev., 

cr 

( 5sin0) 

Probable 
Error  in 
the  Mean, 

a/</n 

Mean 

Probable 
Error  in 
Mean, 

a/Zn 

1 

Mean 

Probable 
Error  in 
Mean, 

cr/Zn 

Applicable  from  May  11,  1961, 

to  June  3, 

1961 

-0.0018 

±0.0009 

±0.0004 

-36.5 

±0.05 

+0.0039 

±0.0012 

±0.0003 

+20.0 

■ 

12 

+0.0003 

±0.0038 

±0.0011 

+0.8 

±2.6 

±0.7 

■ 

52.8 

12 

+0.0154 

±0.0033 

±0.0009 

+32.2 

±2.0 

+8,9 

±0.5 

■ 

20 

12 

-0.0115 

±0.0155 

±0.0045 

-9.1 

±3.5 

±1.0 

■ 

16 

11 

-0.0269 

±0.0226 

±0.0068 

-17.0 

±4.3 

■9 

±1.2 

■ 

Applicable  after  June  12,  1961 

520 

7 

+  0.0007 

±0.0016 

±0.0006 

+15.0 

±12.8 

+4.2 

±3.5 

+0.04 

±0.09 

128 

7 

-0.0026 

±0.0030 

±0.0011 

-13.0 

±5.8 

-3.6 

±1.6 

60 

7 

0.0000 

±0.0063 

±0.0024 

+0.1 

±5.7 

0.0 

±1.6 

52 

7 

-0.0200 

±0.0075 

±0.0028 

-41.1 

±5.8 

±1.6 

20 

6 

-0.0206 

±0.0164 

±0.0067 

-16.3 

±5.3 

±1,5 

16 

7 

+0,0304 

±0.0402 

±0.0152 

+19.3 

±9.6 

+5.4 

±2.7 

*  2 
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is  due  in  part  to  the  limited  precision  of  measurement  used.  TTie  report  part  to  the  limited  precision  of  measurement  used.  Tlie  report 

describes  the  work  done  prior  to  the  second  phase  of  calibration,  that  phase  describes  the  work  done  prior  to  the  second  phase  of  calibration,  that  phase 

leading  to  a  much  higher  precision.  leading  to  a  much  higher  precision. 
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